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Abstract 



We consider the possibility that the cosmological dark matter consists of particles very close in mass to new colored 
particles below the TeV scale. While such a scenario is inherently difficult to directly confirm at colliders, we find 
that indirect dark matter searches may be a powerful alternative. In particular, we show that in this case dark matter 
annihilation to qqg final states can give rise to significant antiproton (but also gamma-ray) fluxes, and compare the 
resulting constraints to bounds from direct searches at LEP, the Tevatron and the LHC. For supersymmetric neutralinos 
degenerate with squarks, e.g., antiprotons can give rise to more stringent constraints for masses below around 45-80 GeV. 

(N 



(N 



£1 
i 

Oh 



(N 
> 
oo 
in 

in 



X 



1. Introduction 

Collider experiments have pushed the scale for pos- 
sible new physics beyond the standard model (BSM) to 
ever higher energies in recent years. The CERN LHC, 
in particular, is now running at a center of mass energy 
of 7 TeV and the non-observation of any clear BSM sig- 
nal, so far, has allowed to place impressive limits p] [2]; 
after an integrated luminosity of 1.04 fb _1 , scalar quarks 
(squarks) in supersymmetric extensions of the standard 
model, e.g., have already been excluded for masses below 
~ 900 GeV OSj. One has to keep in mind, however, that 
these results rest on the assumption that the new colored 
states quickly decay into the lightest neutralino, which is 
assumed to be massless for the sake of the analysis, and 
emit high-energy QCD jets in this process. These limits 
thus do not apply if the only accessible lighter states are 
very close in mass |6j [7l [8] , in which case the resulting 
jets and missing transverse energy would be too soft to 
pass the signal selection criteria. Scenarios with degen- 
erate particle spectra are thus generically very difficult to 
probe at hadron colliders; the situation is considerably bet- 
ter for electron-positron colliders like LEP, but even in this 
case the limits from direct searches for colored states do 
not apply to highly degenerate spectra [HI OH HU HI H3] • 
Here, we investigate whether such a situation could be 
probed by indirect dark matter (DM) searches, assuming 
that the lightest BSM state is a weakly interacting mas- 
sive particle (WIMP) that makes up the cosmological DM. 
The self-annihilation of WIMP pairs could then leave an 
imprint in the spectrum of cosmic rays (see, e.g., Ref. |14| 
for a review on particle DM and indirect searches) . As we 
will demonstrate, WIMP annihilation into qqg final states, 
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a channel so far hardly explored in indirect DM searches, 
may give rise to the most stringent constraints in degener- 
ate scenarios (mostly through antiproton, but also through 
gamma-ray production). In fact, this channel may help to 
fill remaining loop-holes for the existence of new colored 
states below masses of around 100 GeV that are left from 
direct collider searches for such particles. 

The rest of this Letter is organized as follows: In Sec- 
tion 2 we discuss the importance and main characteristics 
of DM annihilation into qqg final states, focussing on the 
case of neutralino DM; in Section 3 we present the energy 
spectra of antiprotons that result from qqg hadronization 
and discuss antiproton propagation through our Galaxy; 
Section 4 contains the limits on the annihilation cross- 
section that can be derived from cosmic-ray observations 
and Section 5 confronts these limits with direct collider 
searches for charged colored particles. Finally, in Section 
6 and 7, respectively, we present a discussion of our results 
and conclude. 

2. Dark matter annihilation into gluons 

Sizable cosmic-ray fluxes from galactic DM can only 
be expected from s-wave annihilation processes; this is be- 
cause the p-wave rate scales like v 2 and is heavily sup- 
pressed for the typically expected DM velocities of the 
order of v ~ 10 -3 . However, if the annihilating DM pair 
is in a J = state — as is necessarily the case for scalar or 
Majorana DM particles x — the tree-level s-wave annihila- 
tion rate into fermion pairs // is helicity suppressed and 
scales like av oc nij/m^. In this situation, the existence of 
an additional vector boson in the final state may lift the 
helicity suppression and lead to a radiative 'correction' of 
the order of (oj/tt) m^/ntf, which in practice can corre- 
spond to an annihilation rate that is enhanced by several 
orders of magnitude [TSl HE] ■ 
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Figure 1: Feynman diagrams contributing to \X ~> 119 m th° limit 
of vanishing quark mass and neutralino velocity. 



We will in the following focus on the supersymmetric 
neutralino, but note that we expect qualitatively very sim- 
ilar results for other DM candidates where the annihilation 
rate into fermions is helicity suppressed. The annihilation 
of neutralinos into qq two-body final states, in particular, is 
typically strongly suppressed (at least for mt, <C m x < m t 
or m x 3> mt), so that the dominant contribution to the 
total annihilation cross section into quarks is given by the 
process XX ~^ QQ9- in the m q — > limit, only the squark- 
exchange i-channel diagrams shown in Fig. [T] contribute 
and we find 




where (i RlL = m 2 qRL /m x and g R P L {g L Pn) is the coupling 
between neutralino, quark and right-handed (left-handed) 
squark; Li 2 (z) = Y^kLi z /k 2 ^ s * ne dilogarithm. This 
result is consistent with the differential cross section pre- 
sented in Ref. [T7] ; for comparison, it can also be obtained 
by the corresponding expressions for photon final states 
997 PS] by simply replacing Q 2 a cm — > (4/3) a s , where Q 
is the electric charge of q^\ 

If the neutralino is a pure Bino with mj, <C m^ < m t , 
the total annihilation cross section is well approximated by 
Eq. (JlJ as long as the mass-splitting /i remains small. This 
is thus the minimal neutralino annihilation cross section 
we can expect; in the exactly degenerate case, fin = hl = 
1, it is largest and becomes 
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For other neutralino compositions, e.g. if the neutralino is 
a Wino, one can find higher total cross sections — mainly 
because annihilation into W + W~ and ZZ final states be- 
comes effective (but also because the neutralino-(s)quark 
couplings in Eq. ([I]) can be larger). However, these sce- 
narios tend to be more constrained than Bino neutralinos, 
and we will here concentrate on the latter. 
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Figure 2: Differential number of antiprotons, per annihilation, from 
various channels and as function of the antiproton kinetic energy T. 
Red (black) curves show the result from three-body (two-body) final 
states containing up-quarks (solid lines) and bottom-quarks (dashed 
lines), respectively. For thick lines, the DM mass is set to 200 GeV, 
for thin lines to 20 GeV. We adopted a mass splitting of fj, = 1.2, but 
the spectra are practically independent of fi in the range fi = 1 . . . 4. 



3. Antiproton production and propagation 

The fragmentation of quarks and gluons into color sin- 
glets leads to the production of a sizable amount of an- 
tiprotons. We derived the antiproton energy spectrum 
dNp/dT (with T denoting kinetic energy) from DM an- 
nihilation using Monte Carlo methods. To this end, we 
simulated the distribution of hard partons from XX ~^ QQ9i 
following our above analytical results, and used the event 
generator Pythia 6.4.19 [20 to perform the subsequent par- 
ton showering, fragmentation and particle decay. 

Our results for the antiproton energy spectra are shown 
in Fig. [2] for different DM annihilation channels and [i = 
1.2 (assuming that /i = /i£ = fiR); for other values of \i 
we obtain very similar results. We find that, in the phe- 
nomenologically relevant low energy region, the additional 
hard gluon in qqg final states leads to an enhancement 
of antiproton production by a factor of up to ~ 2 rela- 
tive to the qq final states commonly considered in indirect 
DM searches. Note that while final state electroweak gauge 
bosons at first sight can lead to a much larger enhancement 
of dNp I dT for low-energy antiprotons (21] , this is only be- 
cause final state radiation of gluons is already included in 
the commonly adopted 'two-body' result obtained from 
Pythia — while that of electroweak gauge bosons is not. 
Fig. [2] also shows that the amount of antiprotons increases 
with decreasing quark mass, as is expected for final state 
radiation (which is dominated by collincar gluons). 

Once produced, antiprotons do not travel along straight 
lines like, e.g., gamma rays, but scatter on randomly dis- 
tributed galactic magnetic field inhomogeneities. Their 
propagation can thus nicely be described in terms of a 
phcnomcnological diffusion model, the free parameters of 



1 Note that there is a typo in Eq. (2) of Ref. [18] — see also Ref. [19] . 
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4. Indirect detection constraints for qqg final states 
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Figure 3: Limits on the DM annihilation cross-section (o~v), as- 
suming a 100% branching ratio into uug final states. Gamma- 
ray limits (dashed line) are taken from Fermi LAT dwarf galaxy 
observations |36j and antiproton limits (solid line) are derived 
from PAMELA data 25 . Dotted lines show the theoretically ex- 
pected values for Bino DM and mass splittings of u = m|/m^ = 
1.01,1.2,2.0 (from top to bottom); the cosmic-ray limits are practi- 
cally independent of fj,. 

which are strongly constrained by other cosmic ray data, 
in particular the boron over carbon ratio B/C [55]. This 
allows to predict the expected astrophysical background, 
consisting of secondary antiprotons mainly produced in 
the collisions of cosmic ray protons with the interstellar 
medium, with a remarkable accuracy |231 124) . The re- 
cent cosmic-ray antiproton flux measurements undertaken 
by the PAMELA experiment over rigidities of 0.35-180 
GV [55] are fully consistent with this background. 

The expected flux of primary antiprotons from DM an- 
nihilation [26], on the other hand, is subject to greater un- 
certainties [55] . The reason for this is that the B/C anal- 
ysis, being restricted to sources in the galactic disk, actu- 
ally leaves large degeneracies in the diffusion parameters — 
in particular between the size of the diffusion constant 
and the thickness L of the diffusion zone perpendicular 
to the galactic plane. While secondary antiprotons are 
not greatly affected by this degeneracy, DM annihilation 
happens predominantly in the galactic halo; therefore, a 
much larger volume of the diffusion zone is probed and 
the antiproton flux is rather sensitive to the adopted value 
of L. Even though the B/C analysis in principle allows 
values as small as L ~ lkpc, values up to L ~ lOkpc are 
preferred when adding radioactive isotopes to the analy- 
sis |27j , with similar conclusions following from the consid- 
eration of gamma rays [28] and cosmic-ray electrons [29] . 
Also radio data are in rather strong conflict with a halo size 
as small as ~ lkpc [3D]. In our analysis, we will therefore 
refer to the 'KRA' model of Ref. [3T] and the B/C best fit 
('MED') model of Ref. [35], both featuring L = 4kpc, as 
well as the 'MAX' model of Ref. [22], featuring L = lOkpc. 



We calculated the antiproton flux from DM annihila- 
tion, as expected at the top of the atmosphere after prop- 
agation, by using DarkSUSY [35] (which implements the 
procedure described in Refs. [3TJ[33]). Solar modulation 
is taken into account by adopting the force-field approxi- 
mation, with a Fisk potential of (j> = 500 MV, but our lim- 
its actually do not strongly depend on this choice because 
they typically do not derive from antiproton energies much 
smaller then 1 GeV. We choose an Einasto profile to model 
the galactic DM halo, with parameters as used in Ref. [3T] ; 
for a Navarro-Frenk- White profile, we get very similar re- 
sults. Furthermore, we assume that the neutralino makes 
up all of the observed dark matter (noting that, due to the 
potentially very efficient coannihilation between squarks 
and neutralinos, this can require a non-thermal produc- 
tion mechanism or a non-standard freeze-out history). 

In order to derive limits on the flux of primary antipro- 
tons from DM annihilation we require that the DM signal 
plus the background of secondary antiprotons does not 
overshoot the antiproton measurements by PAMELA [55] 
by more than 3er in any data point; to be conservative, 
we adopt the minimal astrophysical background flux from 
Ref. [24 . The resulting limits on the annihilation cross- 
section (<rv) are shown in Fig. [3j as function of the DM 
mass m x , for the case of annihilation into uug final states 
and adopting a mass splitting of /x = 1.2; we obtain very 
similar results for other values of /i = 1 ... 4 and for the 
light quark flavors q = d, s, but the limits are up to a fac- 
tor of two weaker in case of the heavier quarks q — c,b. 
We note that these limits are in rather good agreement 
with previous results, like in Ref. [3T], for two-body final 
states. 

Besides antiprotons, the hadronization of quarks and 
gluons also generates a large amount of gamma rays, which 
mostly stem from 7r°-decay. Using again Pythia, we find 
reasonable agreement between the gamma-ray energy spec- 
trum from qqg and the commonly adopted bb spectrum 
(at the peak of the spectra at E 1 ss m x /20 the fluxes 
agree within 20% for q = d,u, s, c, 6). Currently, the Fermi 
Large Area Telescope (LAT) is measuring the gamma-ray 
sky with unprecedented precision and strong limits on the 
DM annihilation cross section were derived from observa- 
tions of dwarf galaxies [311 022 [36] , galaxy clusters [37[ [38] 
and the isotropic diffuse gamma-ray backround 39 . In 
Fig. [3] we show the corresponding limits on the annihila- 
tion cross-section to qqg as derived from the Fermi LAT 
data by a combined analysis of ten dwarf spheroidal galax- 
ies |36] , making use of the above-mentioned similarity be- 
tween the gamma-ray spectra from qqg and bb final states. 
It is reassuring that these limits are almost as strong as 
our antiproton bounds, since limits derived from antipro- 
ton and gamma-ray observations are in principle subject 
to very different astrophysical uncertainties. 

For comparison, we also indicate in Fig. [3] the value of 
the annihilation cross-section for XX ~* uug in the case of 
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Figure 4: Constraints on squarks degenerate with a Bino LSP, as 
function of the Bino mass and the mass splitting. The colored areas 
are excluded by different collider searches as described in the text 
and the areas below the black lines are constrained by cosmic-ray 
antiproton observations (for the 'MED' propagation model). Dotted 
lines correspond to light right-handed up-type squarks, dashed lines 
assume degeneracy between all d, u, c and s-type squarks. The red 
lines show how the limits strengthen in case of the 'MAX' propaga- 
tion model. 



pure Bino DM, as function of the neutralino mass m x and 
the mass splittings \x. Depending on the value of \i, we 
find that DM masses up to ~ 50 GeV can be excluded by 
antiproton and gamma-ray limits; this will be confronted 
with collider limits in the next section. 

5. Collider constraints on new colored particles 

In general it is difficult to directly constrain the DM 
mass with collider experiments. On the other hand, light 
charged and colored particles can be searched for with 
many experiments, and in the recent past the LHC has put 
strong lower limits on squark masses, reaching up to ~ 900 
GeV [3 . However, these constraints are highly model de- 
pendent and, in the case of supersymmetry, depend cru- 
cially on the supersymmetric particle mass spectrum: very 
heavy gluinos, e.g., would suppress the squark production 
in hadron colliders. Additionally, in the parameter re- 
gion in which the mass difference between squarks and 
the lightest superpartner (which is the DM candidate) is 
very small, the constraints from collider experiments are 
evaded since the energy of QCD jets and missing transverse 
momentum that stem from q — > q\ are small, and these 
events therefore cannot pass triggers and cuts. In general, 
for light squarks with masses less than ~ 100 GeV, the 
constrains are very severe because such light particles are 
produced at LEP experiments. Even if the mass splitting 
is relatively small, ~ 5-10 GeV, a large parameter region 
of the squark mass in which it can be pair produced is 
constrained [HI EH HU [HI E3] ■ For squarks that are lighter 
than ~ mz/2, the Z boson decay width generally gives 
the strongest constraints |40) . 

Light squarks are also indirectly constrained by exper- 
iments due to the contribution of loop diagrams; con- 



tributions of these particles have e.g. a large effect on 
electroweak precision measurements and flavor changing 
neutral currents [31] • However, since these limits on the 
squark mass are generically very model-dependent, we will 
not discuss them in this Letter. 

In Fig.[4]we show the constraints on squarks degenerate 
with a Bino LSP that derive from collider null searches, 
compared with the constraints from cosmic-ray antiproton 
observations. In particular, we show the following limits: 
squark and sbottom mass constraints from the LEP exper- 
iments pTJ [13] , sbottom mass constraints from the Teva- 
tron DO experiment [32] and constraints from the LHC 
ATLAS experiment on squark and sbottom masses [51 [33] 
(for constraints from the CMS experiment, see also [33])- 
Concerning the antiproton constraints, we consider two 
limiting cases: (1) only the right-handed up-type squark 
is light (dotted lines), (2) d,U,c and s-type squarks are 
mass degenerate and light (dashed lines) . The areas below 
the corresponding curves are excluded by antiproton mea- 
surements with PAMELA, where we adopted the 'MED' 
propagation model (for the black curves). 

As expected, these limits are strongest for small mass 
splittings and are hence constraining the parameter space 
from a direction that is complementary to the collider lim- 
its. As can be seen from Fig. [4] in case of light up-type 
squarks and small mass-splitting, we can exclude Bino 
masses up to 45 GeV (the 'KRA' propagation model gives 
limits that are about 10 GeV stronger); under favorable 
astrophysical conditions (the 'MAX' model) these limits 
could strengthen to above 60 GeV. If all d, u, c and s-type 
squarks are mass degenerate, the lower limits on the Bino 
mass increase to around 55-80 GeV, depending on the as- 
trophysical scenario. Hence, we find that cosmic-ray obser- 
vations, in particular the observations of antiprotons, can 
be a powerful probe of a parameter region of DM models 
that is generically difficult to access directly with colliders. 

6. Discussion 

While we have demonstrated that favorable astrophysi- 
cal conditions, in particular a large diffusive halo like in the 
'MAX' model, could significantly improve the constraining 
power of indirect DM searches with antiprotons, we would 
like to stress that prospects could be even better when tak- 
ing into account the possibility that the DM distribution 
does not follow a smooth Einasto profile. For an adiabat- 
ically compressed profile as in Refs. [33 [3H], which has 
been argued to result from the gravitational impact of the 
observed galactic distribution of baryons on the DM pro- 
file, e.g., we find a further enhancement of the antiproton 
flux by a factor of about 4-5 at 1 GeV, which corresponds 
to an improvement for the lower bounds on the DM mass 
by about 20-30 GeV. Another rather large enhancement 
could result from the fact that DM is not expected to be 
distributed smoothly, but to cluster in the form of small 
subhalos [37] , The resulting effective 'boost-factor' of the 
antiproton flux could in principle be quite large [24, 48"IB"9"]. 
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though one should note that statistically a much smaller 
value (corresponding to a factor of a few enhancement) 
is expected when extrapolating the results from iV-body 
simulations of gravitational clustering [SU]. Here, we will 
not study such effects in more detail but simply note that 
already a moderate boost-factor of 5 (10) would improve 
the lower limit on m x by roughly 25 (40) GeV. 

Future prospects for indirect DM searches with an- 
tiprotons may be even more promising, given that new 
data with excellent statistics can soon be expected from 
the recently launched AMS-02 experiment [51]. This can 
only lead to better limits on exotic contributions to the 
cosmic-ray antiproton flux, however, if complemented with 
an improved understanding of the astrophysical background 
It was estimated that this could strengthen current bounds 
by about one order of magnitude |31) , which would make it 
possible to probe Bino masses up to 100 GeV and beyond, 
depending on the size of the mass-splitting, the astrophys- 
ical propagation model and the nature of the light squark 
flavors. In addition, as can be seen from Fig.[3j gamma-ray 
observations are also quite constraining for the channels 
considered here; the current limits from Fermi LAT obser- 
vations will strengthen further with a better understanding 
of the Galactic center region and more accumulated data. 

In this Letter we have concentrated on Bino-like neu- 
tralinos, which are probably the most difficult to probe 
due to their relatively weak interactions with fermions and 
their suppressed annihilation signal in cosmic rays. Let 
us stress that the discussed channel qqg can also be rele- 
vant for Wino- or Higgsino-like neutralinos, in which case 
it is however often dominated by other two-body chan- 
nels like, in particular, \X ~> and W + W~ . In these 
cases, masses below a few hundred GeV are difficult to rec- 
oncile with existing gamma-ray or antiproton limits, see 
e.g. Refs. [3TJ[3()]> and we leave a phenomenological explo- 
ration of qqg final states at 0(100 GeV-1 TeV) masses for 
future work. 

Finally, let us mention another important possibility to 
search for degenerate mass spectra. The squark-neutralino 
degeneracy is also expected to yield strong constraints 
from direct dark matter searches because the neutralino- 
nucleus scattering cross section can be resonantly enhanced 
due to the s-channel squark exchange diagram. In fact, 
the XENON100 experiment can have sufficient sensitivity 
for DM masses around 50 GeV [55]. A detailed analy- 
sis of direct detection constraints for any possible squark- 
neutralino degeneracy, however, is beyond the scope of this 
letter and left for future work (but see Ref. [53]). 

7. Conclusions 

In this Letter, we have described a DM annihilation 
channel, and discussed its phenomenological implications 
in some detail, that so far has hardly been explored in 
the context of indirect dark matter searches. If DM is a 
Majorana fermion Xi or a scalar, direct annihilation via 
XX ~~ ^ 99 is strongly helicity suppressed. In this case, the 



overall annihilation rate today can be dominated by inter- 
nal gluon Bremsstrahlung, \X ~^ 995) an d can be boosted 
by several orders of magnitude with respect to the tree- 
level result; while this enhancement is much smaller in the 
early universe, during freeze-out, it can also have a non- 
negligible effect on the relic density JTHl HZJ [S3] . The pre- 
cise branching ratios depend on the difference between the 
DM particle's mass m x and the squark mass rriq; the pro- 
cess XX ~~ ^ 995 becomes strongest in the mass degenerate 
case, i.e. for m,q/m x ~ 1. 

We studied the spectra of cosmic rays generated by 
this channel, and found that the energy spectrum of an- 
tiprotons coming from qqg is generically enhanced with 
respect to the spectrum from qq final states (Fig. [2]). The 
corresponding limits on the annihilation cross-section for 
XX 995 are hence somewhat stronger than limits on the 
traditional bb or W + W~ final states (Fig. From this 
figure we can see, as noted before [3T] ISSTTSl)! [57], that 
antiprotons are most constraining for light DM masses, 
which in particular seems to disfavor the simplest DM in- 
terpretations of recent results from direct DM detection 
experiments [58 ] 159 ] 160]. 

Even more interestingly, the region of small mass split- 
tings between squarks and the LSP (where the annihilation 
fluxes are largest) is precisely the region that is generi- 
cally difficult to probe directly with collider searches for 
squarks; this is because the energy stored in hadronic jets 
coming from the decay q — > \q becomes too small to pass 
the trigger selection in this limit. In this sense, indirect 
DM searches are complementary to collider searches and 
constrain a part of the parameter space that is not directly 
accessible by either LEP, Tevatron or LHC (Fig.|i|(^] De- 
pending on the cosmic-ray propagation model, these limits 
can exclude Bino masses up to 45-80 GeV in the degen- 
erate case. Ongoing AMS-02 measurements are expected 
to strengthen the limits on the annihilation cross-section 
by up to an order of magnitude (see e.g. Ref. |31|). which 
will make it possible to probe Bino masses of 100 GeV and 
beyond. 

Note added. During the final stages of this work, we be- 
came aware of another project analyzing the annihilation 
of DM into qqg final states [5T| . 
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